Stronger genetic differentiation among within-population genetic groups than among populations in Scots pine provides new insights into within-population genetic structuring.

We investigated the presence of spatial genetic groups within forest tree populations and determined if the genetic divergence among these groups is greater than that between populations using Scots pine (Pinus sylvestris) as a model species. We genotyped 890 adult trees of Scots pine in six natural populations in Lithuania at 11 nuclear microsatellite loci. We used a Bayesian clustering approach to identify the within-population genetic groups within each of the six populations. We calculated the differentiation indexes among the genetic groups within each population and among the six populations by ignoring the genetic groups. The Bayesian clustering revealed 2 to 6 distinct genetic groups of varying size as the most likely genetic structures within populations. The genetic differentiation indexes among the genetic groups within populations were nearly tenfold greater (FST = 0.012-0.070) than those between the populations (FST = 0.003). We conclude on the existence of markedly stronger structuring of genetic variation within populations than between populations of Scots pine in large forest tracts of northern Europe. Such genetic structures serve as a contributing factor to large within population genetic diversity in northern conifers. We assume that within population mating in Scots pine is not completely random but rather is stratified into genetic clusters. Our study provides pioneering novel key insights into structuring of genetic variation within populations. Our findings have implications for examining within-population genetic diversity and genetic structure, conservation, and management of genetic resources.

Postglacial phylogeography, admixture, and evolution of red spruce (Picea rubens Sarg.) in Eastern North America.

Climate change is a major evolutionary force that can affect the structure of forest ecosystems worldwide. Red spruce (Picea rubens Sarg.) has recently faced a considerable decline in the Southern Appalachians due to rapid environmental change, which includes historical land use, and atmospheric pollution. In the northern part of its range, red spruce is sympatric with closely related black spruce (Picea mariana (Mill.) B.S.P.), where introgressive hybridization commonly occurs. We investigated range-wide population genetic diversity and structure and inferred postglacial migration patterns and evolution of red spruce using nuclear microsatellites. Moderate genetic diversity and differentiation were observed in red spruce. Genetic distance, maximum likelihood and Bayesian analyses identified two distinct population clusters: southern glacial populations, and the evolutionarily younger northern populations. Approximate Bayesian computation suggests that patterns of admixture are the result of divergence of red spruce and black spruce from a common ancestor and then introgressive hybridization during post-glacial migration. Genetic diversity, effective population size (Ne) and genetic differentiation were higher in the northern than in the southern populations. Our results along with previously available fossil data suggest that Picea rubens and Picea mariana occupied separate southern refugia during the last glaciation. After initial expansion in the early Holocene, these two species faced a period of recession and formed a secondary coastal refugium, where introgressive hybridization occurred, and then both species migrated northward. As a result, various levels of black spruce alleles are present in the sympatric red spruce populations. Allopatric populations of P. rubens and P. mariana have many species-specific alleles and much fewer alleles from common ancestry. The pure southern red spruce populations may become critically endangered under projected climate change conditions as their ecological niche may disappear.

Assembly and Annotation of Red Spruce (Picea rubens) Chloroplast Genome, Identification of Simple Sequence Repeats, and Phylogenetic Analysis in Picea.

We have sequenced the chloroplast genome of red spruce (Picea rubens) for the first time using the single-end, short-reads (44 bp) Illumina sequences, assembled and functionally annotated it, and identified simple sequence repeats (SSRs). The contigs were assembled using SOAPdenovo2 following the retrieval of chloroplast genome sequences using the black spruce (Picea mariana) chloroplast genome as the reference. The assembled genome length was 122,115 bp (gaps included). Comparatively, the P. rubens chloroplast genome reported here may be considered a near-complete draft. Global genome alignment and phylogenetic analysis based on the whole chloroplast genome sequences of Picea rubens and 10 other Picea species revealed high sequence synteny and conservation among 11 Picea species and phylogenetic relationships consistent with their known classical interrelationships and published molecular phylogeny. The P. rubens chloroplast genome sequence showed the highest similarity with that of P. mariana and the lowest with that of P. sitchensis. We have annotated 107 genes including 69 protein-coding genes, 28 tRNAs, 4 rRNAs, few pseudogenes, identified 42 SSRs, and successfully designed primers for 26 SSRs. Mononucleotide A/T repeats were the most common followed by dinucleotide AT repeats. A similar pattern of microsatellite repeats occurrence was found in the chloroplast genomes of 11 Picea species.

Genome-wide identification of PEBP gene family members in potato, their phylogenetic relationships, and expression patterns under heat stress.

BACKGROUND: The phosphatidylethanolamine-binding protein (PEBP) gene family is involved in regulating many plant traits. Genome-wide identification of PEPB members and knowledge of their responses to heat stress may assist genetic improvement of potato (Solanum tuberosum). METHODS AND RESULTS: We identified PEBP gene family members from both the recently-updated, long-reads-based reference genome (DM v6.1) and the previous short-reads-based annotation (PGSC DM v3.4) of the potato reference genome and characterized their heat-induced gene expression using RT-PCR and RNA-Seq. Fifteen PEBP family genes were identified from DM v6.1 and named as StPEBP1 to StPEBP15 based on their locations on 6 chromosomes and were classified into FT, TFL, MFT, and PEBP-like subfamilies. Most of the StPEBP genes were found to have conserved motifs 1 to 5. Tandem or segmental duplications were found between StPEBP genes in seven pairs. Heat stress induced opposite expression patterns of certain FT and TFL members but involving different members in leaves, roots and tubers. CONCLUSION: The long-reads-based genome assembly and annotation provides a better genomic resource for identification of PEBP family genes. Heat stress tends to decrease FT gene activities but increases TFL gene activities, but this opposite expression involves different FT/TFL pairs in leaves, roots, and tubers. This tissue-specific expression pattern of PEBP members may partly explain why different potato organs differ in their sensitivities to heat stress. Our study provides candidate PEBP family genes and relevant information for genetic improvement of heat tolerance in potato and may help understand heat-induced responses in other plants.

Heat-stress-induced sprouting and differential gene expression in growing potato tubers: Comparative transcriptomics with that induced by postharvest sprouting.

Crops face increased risk from heat stress due to climate change. Potato (Solanum tuberosum L.) tubers grown in hot summers often have defects including pre-harvest sprouting ("heat sprouts"). We have used 18 potato cultivars to investigate whether heat stress (HS) conditions alone could cause heat sprouting and dormancy changes in tubers. We also examined transcriptomic responses of potato to HS and whether these responses are like those induced by postharvest sprouting. We demonstrated that HS alone caused heat sprouts and shortened postharvest dormancy period, heat-sprouted tubers became dormant after harvest, and cultivars varied substantially for producing heat spouts but there was no clear association with cultivar maturity earliness. Cultivar Innovator did not show any heat sprouts and still had long dormancy. Dormancy-associated genes (DOG1 and SLP) were downregulated in HS tubers like in postharvest sprouting tubers. We have identified 1201 differentially expressed genes, 14 enriched GO terms and 12 enriched KEGG pathways in response to HS in growing tubers of 'Russet Burbank'. Transcriptomic response of 'Russet Burbank' to HS showed significant similarities to that of postharvest non-HS sprouted tubers. Gibberellin biosynthesis pathway was enriched in heat-stressed tubers and was likely involved in heat sprouting and dormancy release. Heat sprouting and postharvest sprouting shared common candidate genes and had significant similarity in gene expression. Our study has significance for selecting potato cultivars for farming, planning storage and utilization of heat-stressed tubers, identifying sprouting-related genes, understanding heat-stress biology, and breeding heat-tolerant potato cultivars, especially for sustainable potato production under climate change.

Genetic Diversity, Structure and Effective Population Size of Old-Growth vs. Second-Growth Populations of Keystone and Long-Lived Conifer, Eastern White Pine (Pinus strobus): Conservation Value and Climate Adaptation Potential.

Whether old-growth (OG) forests have higher genetic diversity and effective population size, consequently higher conservation value and climate adaptive potential than second-growth (SG) forests, remain an unresolved issue. We have tested the hypothesis that old-growth forest tree populations have higher genetic diversity, effective population size (NE ), climate adaptive potential and conservation value and lower genetic differentiation than second-growth forest tree populations, employing a keystone and long-lived conifer, eastern white pine (EWP; Pinus strobus). Genetic diversity and population structure of old-growth and second-growth populations of eastern white pine (EWP) were examined using microsatellites of the nuclear and chloroplast genomes and single nucleotide polymorphisms (SNPs) in candidate nuclear genes putatively involved in adaptive responses to climate and underlying multilocus genetic architecture of local adaptation to climate in EWP. Old-growth and second-growth EWP populations had statistically similar genetic diversity, inbreeding coefficient and inter-population genetic differentiation based on nuclear microsatellites (nSSRs) and SNPs. However, old-growth populations had significantly higher chloroplast microsatellites (cpSSRs) haploid diversity than second-growth populations. Old-growth EWP populations had significantly higher coalescence-based historical long-term NE than second-growth EWP populations, but the linkage disequilibrium (LD)-based contemporary NE estimates were statistically similar between the old-growth and second-growth EWP populations. Analyses of population genetic structure and inter-population genetic relationships revealed some genetic constitution differences between the old-growth and second-growth EWP populations. Overall, our results suggest that old-growth and second-growth EWP populations have similar genetic resource conservation value. Because old-growth and second-growth EWP populations have similar levels of genetic diversity in genes putatively involved in adaptive responses to climate, old-growth, and second-growth populations may have similar adaptive potential under climate change. Our results could potentially be generalized across most of the boreal and temperate conifer forest trees. Our study contributes to address a long-standing issue, advances research field and knowledge about conservation and ecological and climate adaptation of forest trees.

Development and characterization of Novel EST-based single-copy genic microsatellite DNA markers in white spruce and black spruce.

Due mainly to large genome size and prevalence of repetitive sequences in the nuclear genome of spruce (Picea Mill.), it is very difficult to develop single-copy genomic microsatellite markers. We have developed and characterized 25 polymorphic, single-copy genic microsatellites from white spruce (Picea glauca (Moench) Voss) EST sequences and determined their informativeness in white spruce and black spruce (Picea mariana (Mill.) B.S.P.) and inheritance in black spruce. White spruce EST sequences from NCBI dbEST were searched for the presence of microsatellite repeats. Forty-seven sequences containing dinucleotide, trinucleotide, tetranucleotide and compound repeats were selected to develop primers. Twenty-five of the designed primer pairs yielded scorable amplicons, with single-locus patterns, and were characterized in 20 individuals each of white spruce and black spruce. All 25 microsatellites were polymorphic in white spruce and 24 in black spruce. The number of alleles at a locus ranged from two to 18, with a mean of 8.8 in white spruce, and from one to 17, with a mean of 7.6 in black spruce. The expected heterozygosity/polymorphic information content ranged from 0.10 to 0.92, with a mean of 0.67 in white spruce, and from 0 to 0.93, with a mean of 0.59 in black spruce. Microsatellites with dinucleotide and compound repeats were more informative than those with trinucleotide and tetranucleotide repeats. Eighteen microsatellite markers polymorphic between the parents of a black spruce controlled cross inherited in a single-locus Mendelian fashion. The microsatellite markers developed can be applied for various genetics, genomics, breeding, and conservation studies and applications.

Population genomics for wildlife conservation and management.

Biodiversity is under threat worldwide. Over the past decade, the field of population genomics has developed across nonmodel organisms, and the results of this research have begun to be applied in conservation and management of wildlife species. Genomics tools can provide precise estimates of basic features of wildlife populations, such as effective population size, inbreeding, demographic history and population structure, that are critical for conservation efforts. Moreover, population genomics studies can identify particular genetic loci and variants responsible for inbreeding depression or adaptation to changing environments, allowing for conservation efforts to estimate the capacity of populations to evolve and adapt in response to environmental change and to manage for adaptive variation. While connections from basic research to applied wildlife conservation have been slow to develop, these connections are increasingly strengthening. Here we review the primary areas in which population genomics approaches can be applied to wildlife conservation and management, highlight examples of how they have been used, and provide recommendations for building on the progress that has been made in this field.

Single-Locus versus Multilocus Patterns of Local Adaptation to Climate in Eastern White Pine (Pinus strobus, Pinaceae).

Natural plant populations are often adapted to their local climate and environmental conditions, and populations of forest trees offer some of the best examples of this pattern. However, little empirical work has focused on the relative contribution of single-locus versus multilocus effects to the genetic architecture of local adaptation in plants/forest trees. Here, we employ eastern white pine (Pinus strobus) to test the hypothesis that it is the inter-genic effects that primarily drive climate-induced local adaptation. The genetic structure of 29 range-wide natural populations of eastern white pine was determined in relation to local climatic factors using both a reference set of SSR markers, and SNPs located in candidate genes putatively involved in adaptive response to climate. Comparisons were made between marker sets using standard single-locus outlier analysis, single-locus and multilocus environment association analyses and a novel implementation of Population Graphs. Magnitudes of population structure were similar between the two marker sets. Outlier loci consistent with diversifying selection were rare for both SNPs and SSRs. However, genetic distances based on the multilocus among population covariances (cGD) were significantly more correlated to climate, even after correcting for spatial effects, for SNPs as compared to SSRs. Coalescent simulations confirmed that the differences in mutation rates between SSRs and SNPs did not affect the topologies of the Population Graphs, and hence values of cGD and their correlations with associated climate variables. We conclude that the multilocus covariances among populations primarily reflect adaptation to local climate and environment in eastern white pine. This result highlights the complexity of the genetic architecture of adaptive traits, as well as the need to consider multilocus effects in studies of local adaptation.

Post-glacial phylogeography and evolution of a wide-ranging highly-exploited keystone forest tree, eastern white pine (Pinus strobus) in North America: single refugium, multiple routes.

BACKGROUND: Knowledge of the historical distribution and postglacial phylogeography and evolution of a species is important to better understand its current distribution and population structure and potential fate in the future, especially under climate change conditions, and conservation of its genetic resources. We have addressed this issue in a wide-ranging and heavily exploited keystone forest tree species of eastern North America, eastern white pine (Pinus strobus). We examined the range-wide population genetic structure, tested various hypothetical population history and evolutionary scenarios and inferred the location of glacial refugium and post-glacial recolonization routes. Our hypothesis was that eastern white pine survived in a single glacial refugium and expanded through multiple post-glacial recolonization routes. RESULTS: We studied the range-wide genetic diversity and population structure of 33 eastern white pine populations using 12 nuclear and 3 chloroplast microsatellite DNA markers. We used Approximate Bayesian Computation approach to test various evolutionary scenarios. We observed high levels of genetic diversity, and significant genetic differentiation (F ST = 0.104) and population structure among eastern white pine populations across its range. A south to north trend of declining genetic diversity existed, consistent with repeated founder effects during post-glaciation migration northwards. We observed broad consensus from nuclear and chloroplast genetic markers supporting the presence of two main post-glacial recolonization routes that originated from a single southern refugium in the mid-Atlantic plain. One route gave rise to populations at the western margin of the species' range in Minnesota and western Ontario. The second route gave rise to central-eastern populations, which branched into two subgroups: central and eastern. We observed minimal sharing of chloroplast haplotypes between recolonization routes but there was evidence of admixture between the western and west-central populations. CONCLUSIONS: Our study reveals a single southern refugium, two recolonization routes and three genetically distinguishable lineages in eastern white pine that we suggest to be treated as separate Evolutionarily Significant Units. Like many wide-ranging North American species, eastern white pine retains the genetic signatures of post-glacial recolonization and evolution, and its contemporary population genetic structure reflects not just the modern distribution and effects of heavy exploitation but also routes northward from its glacial refugium.

Highly informative single-copy nuclear microsatellite DNA markers developed using an AFLP-SSR approach in black spruce (Picea mariana) and red spruce (P. rubens).

BACKGROUND: Microsatellites or simple sequence repeats (SSRs) are highly informative molecular markers for various biological studies in plants. In spruce (Picea) and other conifers, the development of single-copy polymorphic genomic microsatellite markers is quite difficult, owing primarily to the large genome size and predominance of repetitive DNA sequences throughout the genome. We have developed highly informative single-locus genomic microsatellite markers in black spruce (Picea mariana) and red spruce (Picea rubens) using a simple but efficient method based on a combination of AFLP and microsatellite technologies. PRINCIPAL FINDINGS: A microsatellite-enriched library was constructed from genomic AFLP DNA fragments of black spruce. Sequencing of the 108 putative SSR-containing clones provided 94 unique sequences with microsatellites. Twenty-two of the designed 34 primer pairs yielded scorable amplicons, with single-locus patterns. Fourteen of these microsatellite markers were characterized in 30 black spruce and 30 red spruce individuals drawn from many populations. The number of alleles at a polymorphic locus ranged from 2 to 18, with a mean of 9.3 in black spruce, and from 3 to 15, with a mean of 6.2 alleles in red spruce. The polymorphic information content or expected heterozygosity ranged from 0.340 to 0.909 (mean = 0.67) in black spruce and from 0.161 to 0.851 (mean = 0.62) in red spruce. Ten SSR markers showing inter-parental polymorphism inherited in a single-locus Mendelian mode, with two cases of distorted segregation. Primer pairs for almost all polymorphic SSR loci resolved microsatellites of comparable size in Picea glauca, P. engelmannii, P. sitchensis, and P. abies. SIGNIFICANCE: The AFLP-based microsatellite-enriched library appears to be a rapid, cost-effective approach for isolating and developing single-locus informative genomic microsatellite markers in black spruce. The markers developed should be useful in black spruce, red spruce and other Picea species for various genetics, genomics, breeding, forensics, conservation studies and applications.

Genetic divergence and signatures of natural selection in marginal populations of a keystone, long-lived conifer, Eastern White Pine (Pinus strobus) from Northern Ontario.

Marginal populations are expected to provide the frontiers for adaptation, evolution and range shifts of plant species under the anticipated climate change conditions. Marginal populations are predicted to show genetic divergence from central populations due to their isolation, and divergent natural selection and genetic drift operating therein. Marginal populations are also expected to have lower genetic diversity and effective population size (Ne) and higher genetic differentiation than central populations. We tested these hypotheses using eastern white pine (Pinus strobus) as a model for keystone, long-lived widely-distributed plants. All 614 eastern white pine trees, in a complete census of two populations each of marginal old-growth, central old-growth, and central second-growth, were genotyped at 11 microsatellite loci. The central populations had significantly higher allelic and genotypic diversity, latent genetic potential (LGP) and Ne than the marginal populations. However, heterozygosity and fixation index were similar between them. The marginal populations were genetically diverged from the central populations. Model testing suggested predominant north to south gene flow in the study area with curtailed gene flow to northern marginal populations. Signatures of natural selection were detected at three loci in the marginal populations; two showing divergent selection with directional change in allele frequencies, and one balancing selection. Contrary to the general belief, no significant differences were observed in genetic diversity, differentiation, LGP, and Ne between old-growth and second-growth populations. Our study provides information on the dynamics of migration, genetic drift and selection in central versus marginal populations of a keystone long-lived plant species and has broad evolutionary, conservation and adaptation significance.

Generation, functional annotation and comparative analysis of black spruce (Picea mariana) ESTs: an important conifer genomic resource.

BACKGROUND: EST (expressed sequence tag) sequences and their annotation provide a highly valuable resource for gene discovery, genome sequence annotation, and other genomics studies that can be applied in genetics, breeding and conservation programs for non-model organisms. Conifers are long-lived plants that are ecologically and economically important globally, and have a large genome size. Black spruce (Picea mariana), is a transcontinental species of the North American boreal and temperate forests. However, there are limited transcriptomic and genomic resources for this species. The primary objective of our study was to develop a black spruce transcriptomic resource to facilitate on-going functional genomics projects related to growth and adaptation to climate change. RESULTS: We conducted bidirectional sequencing of cDNA clones from a standard cDNA library constructed from black spruce needle tissues. We obtained 4,594 high quality (2,455 5' end and 2,139 3' end) sequence reads, with an average read-length of 532 bp. Clustering and assembly of ESTs resulted in 2,731 unique sequences, consisting of 2,234 singletons and 497 contigs. Approximately two-thirds (63%) of unique sequences were functionally annotated. Genes involved in 36 molecular functions and 90 biological processes were discovered, including 24 putative transcription factors and 232 genes involved in photosynthesis. Most abundantly expressed transcripts were associated with photosynthesis, growth factors, stress and disease response, and transcription factors. A total of 216 full-length genes were identified. About 18% (493) of the transcripts were novel, representing an important addition to the Genbank EST database (dbEST). Fifty-seven di-, tri-, tetra- and penta-nucleotide simple sequence repeats were identified. CONCLUSIONS: We have developed the first high quality EST resource for black spruce and identified 493 novel transcripts, which may be species-specific related to life history and ecological traits. We have also identified full-length genes and microsatellite-containing ESTs. Based on EST sequence similarities, black spruce showed close evolutionary relationships with congeneric Picea glauca and Picea sitchensis compared to other Pinaceae members and angiosperms. The EST sequences reported here provide an important resource for genome annotation, functional and comparative genomics, molecular breeding, conservation and management studies and applications in black spruce and related conifer species.

Genetic signatures of natural selection in response to air pollution in red spruce (Picea rubens, Pinaceae).

One of the most important drivers of local adaptation for forest trees is climate. Coupled to these patterns, however, are human-induced disturbances through habitat modification and pollution. The confounded effects of climate and disturbance have rarely been investigated with regard to selective pressure on forest trees. Here, we have developed and used a population genetic approach to search for signals of selection within a set of 36 candidate genes chosen for their putative effects on adaptation to climate and human-induced air pollution within five populations of red spruce (Picea rubens Sarg.), distributed across its natural range and air pollution gradient in eastern North America. Specifically, we used FST outlier and environmental correlation analyses to highlight a set of seven single nucleotide polymorphisms (SNPs) that were overly correlated with climate and levels of sulphate pollution after correcting for the confounding effects of population history. Use of three age cohorts within each population allowed the effects of climate and pollution to be separated temporally, as climate-related SNPs (n = 7) showed the strongest signals in the oldest cohort, while pollution-related SNPs (n = 3) showed the strongest signals in the youngest cohorts. These results highlight the usefulness of population genetic scans for the identification of putatively nonneutral evolution within genomes of nonmodel forest tree species, but also highlight the need for the development and application of robust methodologies to deal with the inherent multivariate nature of the genetic and ecological data used in these types of analyses.

Effects of harvesting of increasing intensities on genetic diversity and population structure of white spruce.

Forest harvesting of increasing intensities is expected to have intensifying impacts on the genetic diversity and population structure of postharvest naturally regenerated stands by affecting the magnitude of evolutionary processes, such as genetic drift, gene flow, mating system, and selection. We have tested this hypothesis for the first time by employing widely distributed boreal white spruce (Picea glauca) as a model and controlled, replicated experimental harvesting and regeneration experiment at the EMEND project site (http://www.emendproject.org). We used two approaches. First, genetic diversity and population structure of postharvest natural regeneration after five harvesting treatments (green tree retention of 75%, 50%, 20%, and 10%, and clearcut) were assessed and compared with those of the unharvested control (pristine preharvest old-growth) in two replicates each of conifer-dominated (CD) and mixed-wood (MW) forest, using 10 (six EST (expressed sequence tag) and four genomic) microsatellite markers. Second, genetic diversity and population structure of preharvest old-growth were compared with those of postharvest natural regeneration after five harvesting treatments in the same treatment blocks in one replicate each of CD and MW forests. Contrary to our expectations, genetic diversity, inbreeding levels, and population genetic structure were similar between unharvested control or preharvest old-growth and postharvest natural regeneration after five harvesting treatments, with clearcut showing no negative genetic impacts. The potential effects of genetic drift and inbreeding resulting from harvesting bottlenecks were counterbalanced by predominantly outcrossing mating system and high gene flow from the residual and/or surrounding white spruce. CD and MW forests responded similarly to harvesting of increasing intensities. Simulated data for 10, 50, and 100 microsatellite markers showed the same results as obtained empirically from 10 microsatellite markers. Similar patterns of genetic diversity and population structure were observed for EST and genomic microsatellites. In conclusion, harvesting of increasing intensities did not show any significant negative impact on genetic diversity, population structure, and evolutionary potential of white spruce in CD and MW forests. Our first of its kind of study addresses the broad central forest management question how forest harvesting and regeneration practices can best maintain genetic biodiversity and ecosystem integrity.

Genetic diversity and differentiation of core vs. peripheral populations of eastern white cedar, Thuja occidentalis (Cupressaceae).

PREMISE OF THE STUDY: Geographically peripheral (marginal) populations are expected to have lower genetic diversity and higher genetic differentiation than geographically core (central) populations as a result of supposedly lower effective population size (N(e)) and higher genetic drift, founder effect, fragmentation, and isolation in peripheral than in core populations. Here we address this issue for a long-lived plant species, eastern white cedar (Thuja occidentalis). METHODS: Genetic diversity and population structure of 13 natural populations of eastern white cedar from its Canadian eastern peripheral and core natural ranges in New Brunswick, Nova Scotia, and Prince Edward Island were studied using six nuclear microsatellite DNA markers. KEY RESULTS: The core populations of eastern white cedar had significantly higher allelic diversity (mean A = 8.83, A(r) = 8.13, A(e) = 4.03) and N(e) (428) than the peripheral populations (A = 6.64, A(r) = 6.15, A(e) = 3.12, N(e) = 198). However, expected heterozygosity was similar in the core (H(e) = 0.64) and peripheral (H(e) = 0.60) populations. Genetic differentiation was significantly higher among the peripheral (F(ST) = 0.089) than among the core (F(ST) = 0.032) populations. No genetic differentiation (F(ST)/Φ(RT) = 0.000) was detected between core and peripheral regions. CONCLUSIONS: Peripheral populations have significantly lower N(e) and genetic diversity in terms of allelic diversity (richness) and significantly higher genetic differentiation than the core populations of eastern white cedar in its Canadian eastern range. However, core and peripheral populations have similar levels of expected heterozygosity. Implications for conservation of eastern white cedar genetic resources are discussed.

Higher fine-scale genetic structure in peripheral than in core populations of a long-lived and mixed-mating conifer--eastern white cedar (Thuja occidentalis L.).

BACKGROUND: Fine-scale or spatial genetic structure (SGS) is one of the key genetic characteristics of plant populations. Several evolutionary and ecological processes and population characteristics influence the level of SGS within plant populations. Higher fine-scale genetic structure may be expected in peripheral than core populations of long-lived forest trees, owing to the differences in the magnitude of operating evolutionary and ecological forces such as gene flow, genetic drift, effective population size and founder effects. We addressed this question using eastern white cedar (Thuja occidentalis) as a model species for declining to endangered long-lived tree species with mixed-mating system. RESULTS: We determined the SGS in two core and two peripheral populations of eastern white cedar from its Maritime Canadian eastern range using six nuclear microsatellite DNA markers. Significant SGS ranging from 15 m to 75 m distance classes was observed in the four studied populations. An analysis of combined four populations revealed significant positive SGS up to the 45 m distance class. The mean positive significant SGS observed in the peripheral populations was up to six times (up to 90 m) of that observed in the core populations (15 m). Spatial autocorrelation coefficients and correlograms of single and sub-sets of populations were statistically significant. The extent of within-population SGS was significantly negatively correlated with all genetic diversity parameters. Significant heterogeneity of within-population SGS was observed for 0-15 m and 61-90 m between core and peripheral populations. Average Sp, and gene flow distances were higher in peripheral (Sp = 0.023, σg = 135 m) than in core (Sp = 0.014, σg = 109 m) populations. However, the mean neighborhood size was higher in the core (Nb = 82) than in the peripheral (Nb = 48) populations. CONCLUSION: Eastern white cedar populations have significant fine-scale genetic structure at short distances. Peripheral populations have several-folds higher within-population fine-scale genetic structure than core populations. Anthropogenic disturbances and population fragmentation presumably have significant effects on fine-scale genetic structure in eastern white cedar. Core populations have higher neighborhood size than peripheral populations, whereas gene flow distances are higher in peripheral than in core populations. The results of our study contribute to the knowledge of poorly-understood spatial genetic structure of core versus peripheral populations in plants. As well, the information is of significance for conservation of genetic resources of eastern white cedar and perhaps of other long-lived forest trees with mixed-mating system.

A high-density genetic linkage map of a black spruce (Picea mariana) × red spruce (Picea rubens) interspecific hybrid.

Genetic maps provide an important genomic resource of basic and applied significance. Spruce (Picea) has a very large genome size (between 0.85 × 1010 and 2.4 × 1010 bp; 8.5-24.0 pg/1C, a mean of 17.7 pg/1C ). We have constructed a near-saturated genetic linkage map for an interspecific backcross (BC1) hybrid of black spruce (BS; Picea mariana (Mill.) B.S.P.) and red spruce (RS; Picea rubens Sarg.), using selectively amplified microsatellite polymorphic loci (SAMPL) markers. A total of 2284 SAMPL markers were resolved using 31 SAMPL-MseI selective nucleotide primer combinations. Of these, 1216 SAMPL markers showing Mendelian segregation were mapped, whereas 1068 (46.8%) SAMPL fragments showed segregation distortion at α = 0.05. Maternal, paternal, and consensus maps consistently coalesced into 12 linkage groups, corresponding to the haploid chromosome number (1n = 1x = 12) of 12 in the genus Picea. The maternal BS map consisted of 814 markers distributed over 12 linkage groups, covering 1670 cM, with a mean map distance of 2.1 cM between adjacent markers. The paternal BS × RS map consisted of 773 markers distributed over 12 linkage groups, covering 1563 cM, with a mean map distance of 2.0 cM between adjacent markers. The consensus interspecific hybrid BC1 map consisted of 1216 markers distributed over 12 linkage groups, covering 1865 cM (98% genome coverage), with a mean map distance of 1.5 cM between adjacent markers. The genetic map reported here provides an important genomic resource in Picea, Pinaceae, and conifers.

Near-saturated and complete genetic linkage map of black spruce (Picea mariana).

BACKGROUND: Genetic maps provide an important genomic resource for understanding genome organization and evolution, comparative genomics, mapping genes and quantitative trait loci, and associating genomic segments with phenotypic traits. Spruce (Picea) genomics work is quite challenging, mainly because of extremely large size and highly repetitive nature of its genome, unsequenced and poorly understood genome, and the general lack of advanced-generation pedigrees. Our goal was to construct a high-density genetic linkage map of black spruce (Picea mariana, 2n = 24), which is a predominant, transcontinental species of the North American boreal and temperate forests, with high ecological and economic importance. RESULTS: We have developed a near-saturated and complete genetic linkage map of black spruce using a three-generation outbred pedigree and amplified fragment length polymorphism (AFLP), selectively amplified microsatellite polymorphic loci (SAMPL), expressed sequence tag polymorphism (ESTP), and microsatellite (mostly cDNA based) markers. Maternal, paternal, and consensus genetic linkage maps were constructed. The maternal, paternal, and consensus maps in our study consistently coalesced into 12 linkage groups, corresponding to the haploid chromosome number (1n = 1x = 12) of 12 in the genus Picea. The maternal map had 816 and the paternal map 743 markers distributed over 12 linkage groups each. The consensus map consisted of 1,111 markers distributed over 12 linkage groups, and covered almost the entire (> 97%) black spruce genome. The mapped markers included 809 AFLPs, 255 SAMPL, 42 microsatellites, and 5 ESTPs. Total estimated length of the genetic map was 1,770 cM, with an average of one marker every 1.6 cM. The maternal, paternal and consensus genetic maps aligned almost perfectly. CONCLUSION: We have constructed the first high density to near-saturated genetic linkage map of black spruce, with greater than 97% genome coverage. Also, this is the first genetic map based on a three-generation outbred pedigree in the genus Picea. The genome length in P. mariana is likely to be about 1,800 cM. The genetic maps developed in our study can serve as a reference map for various genomics studies and applications in Picea and Pinaceae.

A simple method for estimating genetic diversity in large populations from finite sample sizes.

BACKGROUND: Sample size is one of the critical factors affecting the accuracy of the estimation of population genetic diversity parameters. Small sample sizes often lead to significant errors in determining the allelic richness, which is one of the most important and commonly used estimators of genetic diversity in populations. Correct estimation of allelic richness in natural populations is challenging since they often do not conform to model assumptions. Here, we introduce a simple and robust approach to estimate the genetic diversity in large natural populations based on the empirical data for finite sample sizes. RESULTS: We developed a non-linear regression model to infer genetic diversity estimates in large natural populations from finite sample sizes. The allelic richness values predicted by our model were in good agreement with those observed in the simulated data sets and the true allelic richness observed in the source populations. The model has been validated using simulated population genetic data sets with different evolutionary scenarios implied in the simulated populations, as well as large microsatellite and allozyme experimental data sets for four conifer species with contrasting patterns of inherent genetic diversity and mating systems. Our model was a better predictor for allelic richness in natural populations than the widely-used Ewens sampling formula, coalescent approach, and rarefaction algorithm. CONCLUSIONS: Our regression model was capable of accurately estimating allelic richness in natural populations regardless of the species and marker system. This regression modeling approach is free from assumptions and can be widely used for population genetic and conservation applications.